Journal of Animal Ecology

O_LISeasonal forcing is a dominant factor shaping host-pathogen interactions and disease dynamics across many wildlife systems, including species impacted by environmentally transmitted parasites. How seasonality in parasite dynamics translates to the host when the infection and disease impacts operate at different timescales, however, remains poorly understood. C_LIO_LIWe investigate how seasonality shapes sarcoptic mange in bare-nosed wombats, Vombatus ursinus, a disease caused by the environmentally transmitted parasitic mite Sarcoptes scabiei, causing a protracted clinical time-course in the host. Using an empirically informed state-based deterministic model we explore how wombat population trajectories are influenced by (i) seasonal constraints to mite survival and (ii) in context of host-pathogen encounter rates, as measured by the ratio of burrows to wombats. C_LIO_LIWe demonstrate three long-term outcomes of wombat-mange: host and parasite extinction, endemic disease, and disease-free. We find seasonal environments narrow the range of host-pathogen encounter rates that support S. scabiei persistence relative to stable environments, and prevalence and population sizes vary more in seasonal compared to stable environments except under moderate host-pathogen encounter rates when seasonal effects are less apparent. We also find that a protracted infectious period is essential for host-parasite coexistence in the wombat-mange system. C_LIO_LIOur seasonal model results are consistent with field observations, such that mange prevalence in natural populations increases during seasons of longer off-host mite survival. Application of these findings suggest management efforts could reduce host population impacts through disease management in seasons with longer off-host parasite survival or reduce the environmental reservoir through disease management in seasons with shorter off-host survival. C_LIO_LIWe provide novel, mechanistic explanations for distinctive population trajectories that arise from a seasonally forced wildlife disease, including climate factors that operate independently on parasites, host demography, and disparate timescales over which seasonality affects parasites and hosts. Broadly, linking seasonality to long-term population dynamics can improve the predictability and management of wildlife diseases, but requires an understanding of how local intrinsic factors interact with seasonal pressures over time. C_LI

O_LISeasonal migration in birds involves a substantial spatial redistribution of avian biodiversity each year and drives seasonal changes in community composition. Migrants experience different combinations of species interactions over space and time, generating regular disassembly and reassembly of bird communities throughout their annual cycles. However, the effects of seasonal migration on phylogenetic community structure remain poorly understood. C_LIO_LIWe assess spatiotemporal variation in phylogenetic community structure of North American passerines to test how seasonal migration restructures the evolutionary relatedness and dominant assembly mechanisms in bird communities throughout the annual cycle. Using distributional projections, we calculated metrics describing the phylogenetic dispersion of passerine communities each week of the year. We then tested the relationship between seasonal turnover in community phylogenetic dispersion and seasonal variation in species richness and proportion of migratory species. C_LIO_LISeasonal migration, by changing spatial patterns of avian diversity, simultaneously drives a complex continental redistribution of phylogenetic community structure. We find evidence of taxonomic scale dependency to our results, wherein throughout North America, the seasonal influx of migrant passerines yields communities that are overall more phylogenetically clustered, yet also exhibit greater phylogenetic overdispersion at smaller taxonomic scales. C_LIO_LISeasonal shifts in phylogenetic dispersion, though complex, track changes in diversity, manifesting as fluctuations in phylogenetic dispersion between northern and southern regions as seasonal migrants move between these regions. Our findings reveal a dynamic continental landscape of phylogenetic community structure directed by the movements of seasonal migrants. C_LI

Space-for-time Cormack-Jolly-Seber (CJS) models have been developed to estimate survival of migrating animals that are imperfectly detected at spatially discrete sampling stations. However, space-for-time CJS models typically ignore survival over time and age-specific probability of movement, missing the diversity of life history strategies. To address these limitations, we extended space-for-time Cormack-Jolly-Seber models to explicitly incorporate time, age, and individual-level covariates to account for diverse life history strategies. We incorporate a sub-model that includes uncertainty in individuals age using a finite-mixture model. The detection and transition probabilities were parameterized using generalized linear models, facilitating flexible model specification and inclusion of spatiotemporal and individual covariates. We apply the model to detections of juvenile steelhead (Oncorhynchus mykiss) from two populations in the Snake River Basin in Idaho, USA, to estimate trends in survival with respect to population, age, time, and the length of individuals. Additionally, we generalized the model so that it can be applied to other systems and populations with different life history strategies and monitoring infrastructure.

Reproduction in species with parental care involves sustaining a brood of offspring through an energetically demanding period, when shifts in resource availability, weather, predation risk, and parental condition can strongly alter offspring survival. The most extreme outcome is complete brood failure (death of all offspring), which is relatively frequent in many bird species and may occur when conditions cross a viability threshold. Although complete brood failure is important for shaping fitness variation and population dynamics, we have limited understanding of how intra- and interspecific density dependence governs these events, or of how factors such as habitat quality and disease burden contribute to them, because deriving this requires fine-scale, individual-level data collected across generations for multiple overlapping species. Using a dataset totalling 38,509 nesting attempts from great tits (Parus major) and blue tits (Cyanistes caeruleus) in Wytham Woods, Oxford, UK, we examined how brood failure is shaped by local conspecific and heterospecific density, habitat structure, and avian malaria infection for a subset. Complete brood failure was frequent (14.75%), mostly involving chick mortality in the nest consistent with starvation, rather than brood removal by predators. Relationships between density and brood failure were strong but species-specific. Specifically, great tit failure risk was higher in neighbourhoods that remained densely populated across years, whereas blue tit failure risk was lower where annual great tit or combined density was high, but not where annual blue tit density itself was high. This suggests that local overall density reflects continuing constraint for great tits, while local annual density may partly track favourable within-year conditions and settlement patterns for blue tits. In great tits, failure was also more common where oak density was low and farther from the closest river (Thames), while habitat associations were weak in blue tits. Malaria infection was spatially heterogeneous and covaried with density and habitat, but infection status did not significantly explain complete brood failure. Together, these results show that complete brood failure is shaped by spatially structured local ecological context, and how density dependence in these events can differ in direction and timescale between sympatric species.

Wildlife diseases cause well-documented and often dramatic reductions in host survival. However, the impact of infectious diseases on host reproduction remains understudied, especially with respect to effects of prior and/or current pathogen exposure on reproductive development. Here we experimentally tested how prior and/or current infection with a common bacterial pathogen, Mycoplasma gallisepticum ( MG), alters reproductive development for female versus male house finches (Haemorhous mexicanus). Finches were inoculated with either MG or sterile media while in wintering condition and subsequently received one of these treatments while in breeding condition. In females, MG exposure had both immediate and carry-over effects on reproduction: controls had higher odds of laying eggs compared to females inoculated with MG in spring only, higher odds than females inoculated in both winter and spring, and higher odds than females given MG in the winter only. Conversely, breeding-condition males inoculated with MG in spring had higher testosterone levels than males receiving only control inoculations, and there were no carryover effects of winter MG inoculation or inoculations during both seasons on testosterone. Sex bias in the reproductive impacts of infectious diseases may have important knock-on effects on the epidemiology and population-regulating effects of pathogens, thereby warranting further study.

ContextThe movement of vertebrate hosts across urbanized landscapes can play a key role in the transmission of direct-contact diseases. Understanding how wildlife hosts move in urban landscapes, and how transmission is affected by their landscape-constrained and disease-altered movements, is imperative for better predicting the spread of disease. ObjectiveWe assess how the movement of red foxes (Vulpes vulpes) according to landcover type, and their infection status, affect the spread of mange (caused by Sarcoptes scabiei) in an urbanized landscape. MethodsWe developed a mange transmission model (MTM) using an agent-based model to compare two movement behaviours of foxes in Scarborough (Ontario, Canada): random and landcover-based. We further assessed the effects of movement on disease transmission by considering the foxs infection status and comparing a range of movement probability scenarios. We quantified the number of effective contact events and the effective reproduction number (Re) according to each scenario. ResultsWe found that both landcover-dependent movement and infection status influenced the spread of mange within fox populations. The number of effective contact events and effective reproduction number Re was greatest when landscape heterogeneity was included in the model and foxes moved through paths of least resistance to movement, and when susceptible and infected foxes had an equal probability of leaving a fragmented habitat patch. ConclusionsOur findings suggest that mange spread may be accelerated along movement corridors in fragmented, heterogenous landscapes. As urban areas expand and remnant habitat within these is further lost and animals are relegated to fewer movement pathways, disease transmission may increase.

In many ecosystems, anthropogenic warming is reshaping thermal regimes, leading to resource quality declines and imposing a dual constraint for ectotherms: elevated metabolic demand coupled with reduced assimilable energy. We tested whether plasticity in gut morphology and gut microbiome can buffer amphibian larvae against these concurrent stressors. Common frog (Rana temporaria) tadpoles were reared at two temperatures (18 vs. 24.5{degrees}C) crossed with three food-quality treatments (low, medium, high). We quantified growth and developmental rates, critical thermal limits (CTmax, CTmin), gut morphology (mass, relative length), and gut bacterial diversity and composition, together with predicted functional pathways. Warming accelerated growth and development and increased CTmax. Food quality increased growth and development, with temperature-dependent effects on developmental rate and CTmax. Gut mass declined at higher temperature and low-quality diets, but relative gut length showed only modest diet effects and no temperature dependence. Bacterial community composition and structure shifted with temperature and food quality. Predicted pathways suggest functional reconfiguration under warming and low food quality, consistent with sustaining energy acquisition and mitigating metabolic and oxidative stress. Together, these results implicate microbiome plasticity, rather than gut morphological plasticity, as a candidate mechanism supporting larval performance and heat-tolerance acclimation under warming and low food quality.

Mortality causes and survival rate estimates from individual follow-ups of GPS-collared animals provide reliable demographic data, but such estimates are rare for red foxes Vulpes vulpes. We used data from 126 GPS-collared red foxes tracked between 2011 and 2019 across a latitudinal gradient in central Sweden and Norway to quantify mortality causes and estimate annual survival probabilities using the Andersen-Gill extension of the Cox proportional hazards model. Hunting was the dominant mortality cause (63%), followed by vehicle collisions (14%), stress or malnutrition (11%), sarcoptic mange (9%), and predation (3%). Annual survival was 0.58 for adults and 0.32 for subadults. Subadults had approximately twice the hazard of adults, males had elevated hazard relative to females, and mortality risk was highest in autumn and early winter. These collar-based estimates are strikingly similar to mark-recovery estimates from the same region nearly five decades earlier, despite an intervening population collapse from sarcoptic mange, subsequent full recovery, and a major decline in harvest pressure. This convergence is consistent with density-dependent regulation at a food-determined and food-limited carrying capacity, suggesting that current harvest levels are not limiting the population. Furthermore, human-caused mortality was approximately four times greater than natural mortality, illustrating that living in close proximity to humans -- which often favours generalist predators such as red foxes -- may also come at a cost of increased mortality risk. Developing validated census methods for red foxes in boreal Scandinavia is identified as a key priority for quantitative population management.

The impacts of increasingly frequent and intense heatwaves on parasitism are an important frontier for understanding disease risk under climate change. These impacts are complex because parasitism arises from multiple interacting host and parasite traits that can vary in thermal sensitivity and among populations adapted to different temperature regimes. Here, we used a lab microcosm experiment to investigate the effects of heatwaves occurring during two different phases of a winter-adapted mosquito host - ciliate parasite interaction, for six pairs of sympatric host and parasite populations sourced from two geographic regions with differing histories of winter heat. We found that because heatwaves allowed mosquito larvae to evade infection, they reduced parasitism and increased survival. An early heatwave during initial parasite attack had stronger effects than a later heatwave occurring after infections had established. We did not find evidence of local adaptation to heatwaves: impacts were consistent regardless of population, and were mechanistically predictable from previously measured thermal performance curves that described lower infection and stronger host defenses at warm constant temperatures. The results suggest that increasingly frequent heatwaves may accelerate geographic shifts in parasitism, and demonstrate how fundamental host - parasite thermal biology links to the impacts of extreme temperature events.

Population dynamics and their component vital rates may be driven by weather, climate teleconnections between sea and air (e.g. ENSO), or biotic interactions. These drivers operate directly or indirectly and on different temporal scales. We used a Bayesian structural equation model to characterize effects among weather, climate, and incidental intraguild predation (IGP) on the butterfly Euphydryas gillettii's vital rate of pre-diapause survival, using an 18 year dataset. IGP was a major determinant of pre-diapause survival, along with direct and indirect effects of weather and spring climate teleconnections. The direction of climate effects was reversed when mediated through IGP. Our analysis illustrates the need for sequential hypotheses to capture the cascading effects of abiotic factors via biotic interactions. Using sequential hypotheses addresses the debate on weather -- climate teleconnection roles by disentangling their contributions from one another. Finally, vital rates must be decomposed to component rates in order to detect their drivers.

Rising temperatures and unpredictable weather patterns are directly linked to emerging infectious diseases that threaten global biodiversity. Thermal variability affects host-parasite interactions, and understanding how animal hosts respond to temperature variability and parasite infection is paramount for conservation and for predicting zoonotic spillover. Hosts employ two strategies to defend against parasites: resistance (inhibiting or limiting infection) and tolerance (limiting the negative effects of infection). This observation raises a fundamental question: How does thermal variability affect the expression and evolution of resistance versus tolerance? Subsequently, how does the evolution of resistance versus tolerance influence parasite load dynamics? Here, we first review why temperature may differentially affect resistance and tolerance mechanisms, thereby altering host selection and eco-evolutionary disease outcomes. Second, to highlight the importance of these interactions, we present a model that illustrates key potential effects of temperature variability on host defense mechanisms. Our model demonstrates that temperature variability alone could drive lower infection prevalence and loads, but it also selects for host tolerance, ultimately leading to higher net prevalence and loads. We also find widely divergent outcomes depending on how temperature impacts defense strategies. These results highlight key areas for future empirical and theoretical work on the interactions among temperature variability, infection load, and host defense evolution.

O_LILocal environmental conditions during the breeding season can limit bird populations. Identifying which variables, when, and how they affect key biological traits such as body condition is crucial for understanding long-term population trends under ongoing climate change. C_LIO_LIWe analysed the relationships between environmental variables and body condition during the breeding season in European robins (Erithacus rubecula), aiming to uncover links between short-term environmental influences and long-term trends in body condition in the context of local climate change. C_LIO_LIUsing data from a robin population monitored between 2007 and 2021, we applied weather sliding-window analyses to identify periods when temperature, soil moisture, and vegetation productivity best predicted individual body condition. For each variable, we identified critical time windows (CTWs) influencing (1) body condition across the season and (2) individual changes within two weeks. Juveniles and adults were analysed separately, with adult males and females distinguished during pre- and post-fledging periods. We also assessed long-term trends in environmental variables and body condition, and examined how body condition was correlated with apparent survival. C_LIO_LIBody condition variation across the season was explained by different environmental variables depending on age, sex, and period. Body condition declined with increasing minimum temperatures in adult males and juveniles, and with low soil moisture in adults of both sexes. We did not identify reliable CTWs explaining short-term within-individual changes in body condition. Across 2007-2021, body condition in adult males during the post-fledging period declined with rising minimum temperatures, while fledging dates advanced. Apparent survival was positively associated with body condition only in juvenile robins. C_LIO_LIOur results reveal multiple seasonal environmental influences that may contribute to short- and long-term declines in body condition in European robins, with effects particularly strong (or most detectable) in adult males. Reduced body condition may have demographic consequences by lowering juvenile survival, although shifts in breeding phenology could mitigate this impact. Overall, these findings highlight how environmental effects on body condition can shape long-term population trends and species vulnerability to climate change. C_LI

Disentangling causation from correlation is the foundation of the scientific method. Yet, growing evidence suggests that much observational ecology research has not correctly made this distinction due to inappropriate statistical modelling and unappreciated time-delays. Here, we apply time-lagged causal inference techniques to assess the drivers of seabird declines, using multi-decadal North Atlantic seabird data across the behaviour, mass, survival, reproduction and population size of two species ecology (Atlantic puffin, Fratercula arctica, and black-legged kittiwake, Rissa tridactyla). We demonstrate that both climate and anthropogenic activity can suppress breeding success and cause population declines. Moreover, population size is specifically impacted by delayed recruitment effects where both species decline after a lag corresponding to their estimated age of first reproduction. These North Atlantic seabirds are therefore at risk from future environmental and anthropogenic changes, as time-delays may result in populations already on an extinction trajectory prior to changes being detectable in their abundance.

Vaccination of free-roaming dogs (Canis lupus familiaris) and cats (Felis catus) remains a major public health challenge. Rapid urbanization forces these species into complex contact zones, where structural failure of pulsed vaccination under high demographic turnover undermines standard One Health interventions. In such cases species-specific intervention cycles are needed to reduce zoonotic disease risk. We integrated field-data within a simulated vaccination campaign to determine how species-specific turnover rates drive the erosion of herd immunity at an Amazonian urban sentinel site (university campus). We monitored free-roaming populations (72 dogs, 75 cats) using a non-invasive photographic mark-resight protocol from 2023 to 2025. We modelled time-to-disappearance using Cox Proportional Hazards and simulated the trajectory of effective vaccination coverage against a 40% herd immunity threshold, distinguishing between loss of vaccinated individuals and recruitment of susceptible individuals. The campus functioned as a high-turnover system, with 72% of dogs and 48% of cats classified as transients. Species significantly predicted persistence (Hazard Ratio = 0.56; 95% CI: 0.33-0.94; p=0.029), with cats exhibiting double the median residency of dogs (432 vs. 193 days). Consequently, the species experienced divergent epidemiological failure modes. For dogs, simulated vaccination coverage collapsed below a 40% herd immunity threshold in 160 days, driven by rapid immunity attrition (the loss of vaccinated individuals). Although cats persisted longer, their effective coverage was eroded by immunity dilution due to recruitment of naive juveniles, creating a 33% gap between cohort survival and population-level immunity by day 365. Annual vaccination campaigns are likely insufficient in this high-turnover urban dog population. Effective One Health zoonotic control strategies must transition from static abundance-based targets to dynamic, species-specific and turnover-adjusted intervention schedules.

Body size is a fundamental determinant of ecological performance, and variation in length- weight relationships can provide insight into how populations respond to local environmental conditions. Annual killifishes inhabit highly seasonal wetlands characterized by substantial environmental heterogeneity, yet population-level variation in growth allometry remains poorly understood. We investigated growth allometry in 18 populations of two annual killifish species, Nothobranchius kirki and N. wattersi, distributed across ephemeral habitats in Malawi. Length- weight relationships were analyzed separately for each species using linear mixed-effects models, with locality incorporated as a random effect. Environmental variation among localities was summarized using principal component analysis of water quality variables, and population- specific allometric coefficients were subsequently related to environmental gradients. Significant population-level variation in growth trajectories was detected in both species, indicating divergence in allometric scaling among localities. Divergence was more pronounced in N. wattersi, which exhibited a broader range of population-specific allometric coefficients than N. kirki. Populations also differed significantly in relative body condition after accounting for body length and sex. Environmental gradients explained a significant proportion of variation in allometric slopes in N. kirki, whereas no significant relationship was detected in N. wattersi. Thus, the species exhibiting weaker allometric divergence showed stronger environmental associations, while the species exhibiting greater divergence showed little correspondence with measured environmental variables. These results demonstrate substantial spatial heterogeneity in growth allometry among populations of annual killifishes inhabiting seasonal wetlands. Furthermore, the contrasting environmental associations observed between species suggest that population divergence in growth trajectories may arise through different ecological and evolutionary processes, even among closely related taxa occupying similar habitats.

The black-legged tick (Ixodes scapularis), a key vector of Lyme disease, anaplasmosis, and babesiosis, exhibits regionally distinct patterns of seasonal activity driven by climate. Consequently, the relative timing of larval and nymphal activity varies across geographic locations, influencing pathogen transmission dynamics. Early-emerging nymphs may increase pathogen transmission, whereas early-emerging larvae may reduce transmission. In addition, synchrony between the two life-stages facilitates co-feeding transmission, which contributes to pathogen maintenance and coinfection risk. Temperature is thought to be an important driver of tick phenology, but existing mechanistic models that incorporate temperature fail to accurately capture the timing of larval and nymphal tick activity. To address this limitation, we developed a mechanistic model that includes two additional factors: humidity-dependent questing and low rates of overwinter development. To assess the value of these factors for explaining real-world patterns, we fitted alternative models to tick collection data from the National Ecological Observatory Network. In doing so, we found that explicitly incorporating humidity is necessary to reproduce observed tick phenology, with larval ticks being especially sensitive to relative humidity compared to other life stages. In addition, we found that accounting for humidity had a larger effect at Mid-Atlantic sites than at Northeastern sites, underscoring the importance of region-specific interactions between temperature and humidity in shaping I. scapularis phenology. By more accurately capturing tick seasonality compared to existing mechanistic models, our model illustrates the importance of accounting for factors beyond temperature for investigating how climate variability influences seasonal tick activity and pathogen transmission.

Arctic and Boreal raptor communities will continue to be affected by borealization and other climate change related processes, providing a challenge for ecologists predicting future sates. However, by using community assembly theory and species traits, future communities may be predictable. In this study, we analyzed variation in reproduction traits as a consequence of diet specialization for 29 raptors, 2 skuas and 3 corvids. We assessed and implemented foraging traits for specialists and generalists into predator-prey models from which successful invasion conditions were derived. Specialist raptors produced larger clutch sizes, had a higher proportion of fledged per clutch and also expressed more variation compared to generalist raptors. These results suggest a relationship between diet specialization and reproductive traits which was also observed within phylogenetic orders. Specialist owls (Strigiformes) produced higher clutch sizes with a larger clutch range compared to generalist owls. The same pattern was observed for falcons (Falconiformes). No clear difference in reproduction was observed for specialist and generalist hawks, kites and eagles (Accipitriformes). Corvids expressed clutch sizes similar to that of specialist raptors while having the lowest proportion of fledged per clutch. Differences in foraging traits between specialists and generalists could be distinguished using functional response curves. A predator-prey model parameterized with foraging trait data showed that a generalist can coexist with a resident specialist if it has access to prey unavailable to the resident specialist. Otherwise, the native specialist outcompetes the invading generalist due to foraging efficiency. The combined empirical and theoretical findings in this study show how diet specialization affects both reproduction and the potential invasion success of raptors.

Age-dependent survival is central to understanding population dynamics and life-history evolution. We analysed carcass weight and age-at-harvest data from 6022 red foxes (Vulpes vulpes) collected across Sweden between 1967 and 1971 to evaluate latitudinal effects on body mass and age-dependent survival. Carcass weights decreased from south to north in both adults and sub-adults, contrary to Bergmann's rule, with southern foxes weighing approximately 1.27 times more than northern foxes. The latitudinal weight gradient exceeded the sex difference in both age classes, and no sex x region interaction was detected. The decrease in weight with latitude is consistent with reduced prey availability and harsher winter conditions in the north, which limit growth and body size during development. Using a Bayesian age-at-harvest model with region-specific population growth rates (lambda), we estimated age-dependent survival probabilities for four latitudinal regions and both sexes. Despite the strong latitudinal gradient in weight, survival did not show a corresponding pattern - regional differences were uncertain, with all credible intervals spanning zero. Regional population growth rates were consistent with slight decline in the north and near-stability in the south-central region, which suggests that body condition and population dynamics are coupled at the regional scale despite no survival gradient. The decoupling of body condition and survival across regions suggests that mortality patterns are similar across the latitudinal gradient. We discuss these patterns in terms of latitudinal productivity gradients, prey availability, and life-history trade-offs in a widely distributed carnivore.

BACKGROUND: Varroa destructor is the major ectoparasite of honey bees and a vector of viral pathogens. Because pathogen transmission and exposure to contact-active acaricides depend on mite host contacts, understanding the factors governing host residence time is important for both disease epidemiology and pest management. We quantified host residence time under varying bee densities and host-type compositions. RESULTS: Mean residence time was 9.48 h across 312 host-residence events. Mites remained on individual hosts for only 2.36 h on Day 1 but approximately 11-14 h from Day 2 onward. A generalized linear mixed model showed a strong positive effect of day on residence time ({beta} = 0.341, SE = 0.044, P < 0.001), corresponding to an approximately 41% increase in residence time per day. Excluding Day 1 eliminated this effect (P = 0.16), indicating that the temporal pattern was driven primarily by the initial exposure period. Reconstructing Day 1 observations to an 8-hour schedule confirmed that this pattern was not an artifact of observation frequency. Neither host type nor bee density affected residence time, and mite occupancy of nurse bees matched host availability. CONCLUSION: Host residence time was governed primarily by initial exposure rather than host identity or moderate crowding. The results identify a previously undescribed exploratory phase immediately after mites enter a novel adult-bee population. Because shorter residence times imply more frequent host switching, these findings improve our understanding of pathogen transmission dynamics and may help explain variation in the performance of contact-based Varroa control strategies.

Carotenoid- and melanin-based plumage coloration traits are key signals in avian communication and sexual selection as they are often thought to provide "honest" information about individual condition and fitness. These traits arise through distinct but interconnected physiological and genetic pathways. Recent work suggests that there may be a link between host-associated gut microbiota and the functional pathways leading to pigment-based plumage coloration, but this remains largely unexplored in wild populations. To address this gap, we tested whether variation in plumage coloration, as well as breeding condition, is associated with gut microbiome variation in wild populations of male Common Yellowthroats (Parulidae: Geothlypis trichas). We quantified multiple plumage coloration traits and characterized gut microbiome bacterial diversity using 16S rRNA metabarcoding. Through a comprehensive modeling framework, we found that individuals with brighter, more orange-tinted breast feathers and smaller cloacal protuberances (a proxy for breeding condition) exhibited higher gut microbiome diversity. At the taxonomic level, Methylobacterium-Methylorubrum, a carotenoid-producing bacteria, showed strong associations with multiple plumage traits, including mask area, breast feather hue, and saturation. Our results demonstrate that gut microbiome diversity is associated with variation in carotenoid-based coloration traits and breeding condition in Common Yellowthroats. More broadly, these results highlight the potential for host-microbiome interactions to shape phenotypic variation through physiological pathways in wild animal populations.